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Abstract 

We derive the two-loop QCD helicity amplitudes for the processes £q — > £qg [Iq — » iqg) and Ig — > Iqq, 
which are the partonic reactions yielding (2 + l)-jet final states in deep inelastic lepton nucleon scattering. 
The amplitudes are obtained by analytic continuation of the known helicity amplitudes for e + e~ — » qqg. 
We separate the infrared divergent and finite parts of the amplitudes using Catani's infrared factorization 
formula. The analytic results for the finite parts of the amplitudes are expressed in terms of one- and 
two-dimensional harmonic polylogarithms. To evaluate these functions numerically, we list in detail the 
non-trivial (and kinematic region dependent) variable transformations one needs to perform. 



1 Introduction 



The study of the hadronic final state in deep inelastic lcpton-nucleon scattering (D1S) allows the determination 
of aspects of the nucleon structure which are not accessible in inclusive scattering as well as testing predictions 
of quantum chromodynamics (QCD) to a high accuracy. In particular, the two-plus-one-jet production [1] 
rate in DIS can be used both to probe the gluon distribution and to enable a precise determination of the 
strong coupling constant a s . 

The measurement [4, 5] of these jet observables at the electron-proton collider HERA, with a much larger 
kinematical range than at earlier fixed target experiments, allows precision studies, and the inclusion of 
next-to-leading order (NLO) corrections [2, 3] has become mandatory. However, the present data already 
highlight the limitations of the NLO description: the error on the extraction of a s from HERA (2+1) jet 
data is dominated [6] not by the statistical uncertainty on the data, but by the uncertainty inherent to the 
NLO calculation, as estimated by varying renormalization and factorization scales. Given that the statistical 
precision of the data will further improve once all data from the HERA-II run are analysed, the theoretical 
description requires the inclusion of the next-to-next-to- leading order (NNLO), i.e. O(a^), corrections. 

The dijet production cross section in deep inelastic scattering at HERA is very sensitive to the gluon 
distribution in the proton. However, present determinations of parton distribution functions [7] at NNLO 
accuracy [8] can only include data sets for observables where NNLO corrections are known [9] . Consequently, 
the precise HERA deep inelastic dijet data could be used in global NNLO determinations only once the NNLO 
QCD corrections to this observable are computed. 

In terms of matrix elements, the calculation of these corrections requires the computation of three contri- 
butions: the tree level j*i — > 4 partons amplitudes, the one-loop corrections to the j*i — ► 3 partons amplitudes 
and the two-loop corrections to the j*i — > 2 partons amplitudes, where i = q,g. All of these amplitudes can 
be obtained by analytic continuation from the respective amplitudes contributing to e + e~ — > 3 jets at NNLO: 
the tree level 7* — > 5 partons [10], the one- loop corrections to the 7* — ► 4 partons amplitudes [11], and the 
two-loop corrections to the 7* — > 3 partons amplitudes [12,13]. The analytic continuation of the tree level 
and one-loop amplitudes is rather straightforward, it has already been carried out in the context of the NLO 
corrections to the (3 + l)-jet rate in DIS [14]. The analytic continuation of the two-loop amplitudes is more 
involved and is the topic of the present paper. 

Corrections to the deep inelastic (2 + l)-jet cross section, which are the phenomenologically most relevant 
applications, require only the computation of the helicity averaged squared matrix element. However, the 
hclicity amplitudes derived here can be applied to compute angular correlations between the outgoing lepton 
direction and the orientation of the (2 + l)-jet system [15]. They also allow the calculation of jet cross 
sections in polarized lepton-nucleon collisions [16] which can constrain the polarized parton distributions in 
the nucleon. 

Even though the matrix elements for all NNLO partonic subprocesses contributing to (2 + l)-jet final states 
in DIS are now available, additional work is still needed before quantitative predictions can be made. The two- 
loop amplitudes computed here must be combined with the one-loop corrections to — > 3 partons, where 
one of the partons becomes collinear or soft, as well as tree-level processes j*i — ► 4 partons with two soft or 
collinear partons in a way that explicitly allows all of the infrared singularities to cancel one another. This task 
is usually accomplished by a subtraction method, which exploits the universal factorization properties of QCD 
amplitudes in infrared soft or collinear limits to construct real radiation subtraction terms. These subtraction 
terms approximate the full real radiation matrix elements in all unresolved limits, but are sufficiently simple 
to be integrated analytically. Consequently, they can be combined with the virtual loop corrections to the 
matrix elements to yield infrared finite results for suitably defined observables. Several general subtraction 
schemes have been proposed at NLO [17]. 

For observables with only final state partons, an NNLO subtraction formalism, antenna subtraction, 
has been established in [18]. The antenna subtraction formalism constructs the subtraction terms from an- 
tenna functions, which are derived systematically from physical matrix elements [19]. This formalism has 
been applied in the computation of NNLO corrections to three-jet production in electron-positron annihila- 
tion [20-22] and related event shapes [23,24]. This formalism can be extended to include parton showers at 
higher orders [25] , thereby offering a process independent matching of fixed-order calculations and logarithmic 
resummations [26,27], which had to be done on a case-by-case basis for individual observables [28] up to now. 
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For processes with inital-state partons, antenna subtraction has been fully worked out only to NLO so far [29], 
work towards an extension to NNLO is in progress. 

Other approaches to perform NNLO calculations of exclusive observables with initial state partons are 
the use of sector decomposition and a subtraction method based on the transverse momentum structure of 
the final state. The sector decomposition algorithm [30] analytically decomposes both phase space and loop 
integrals into their Laurent expansion in dimensional regularization, and performs a subsequent numerical 
computation of the coefficients of this expansion. Using this formalism, NNLO results were obtained for Higgs 
production [31] and vector boson production [32] at hadron colliders. Both reactions were equally computed 
independently using an NNLO subtraction formalism exploiting the specific transverse momentum structure 
of these observables [33] . 

This paper is structured as follows: in Section [2J we describe the kinematical situation relevant to DIS- 
(2 + l)-jet production. The two-loop helicity amplitudes are computed by analytic continuation from the 
7* — > 3 partons amplitudes in Section[3] Finally, SectionQ]contains a discussion of the results and conclusions. 

2 Kinematics 

The partonic subprocesses contributing to the production of two hard jets (besides the current jet) in DIS 
are lepton-quark scattering 

Hps) + q{~P2) — > l(-p») + q(pi) + g(p 3 ) , (2.1) 
(with lepton-antiquark scattering following trivially from charge conjugation) and lepton-gluon scattering 

t{pn) + g{-ps) — > + q(pi) + q{P2) ■ (2.2) 

The momentum of the vector boson mediating the interaction is given by 

-<Z" = -j#-j£. (2-3) 

It is convenient to define the invariants 

S12 = (Pi +P2) 2 , s 13 = (pi +p 3 ) 2 , s 23 = (pa +p 3 ) 2 , (2-4) 

which fulfil 

q 2 = (pi +P2+ P3) 2 = Sl2 + Sis + S 2 3 = S123 , (2.5) 
as well as the dimensionless invariants 

X = S12/S123 , y = Sl3/si23, Z = S 2 3/si23, (2.6) 

which satisfy x + y + z = 1 . 

In e + e~ — > 3 jet (3j) production, q 2 is time-like (hence positive) and all the Sij are also positive, which 
implies that x, y, z all lie in the interval [0; 1], with the above constraint x + y + z = 1. For DIS- (2 + l)j 
production, q 2 is space-like (hence negative) and for lepton-quark scattering (|2.ip the invariants fulfil 

q 2 < , s 23 > , aia < 0, sis < , (2.7) 

or, equivalently, 

x>0, y>0, z<0. (2.8) 
For lepton-gluon scattering (|2.2|) . one finds 

q 2 < , si2 > , si 3 < 0, s 23 < , (2.9) 

respectively, 

x<0, y > , z>0 . (2.10) 

It is worth noting that the kinematical regions for these two processes are further subdivided by cuts in 
x = 1 and y = 1 (lepton-quark scattering) or y — 1 and z — 1 (lepton-gluon scattering), as was first pointed 
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out by Graudenz in the context of the calculation of the NLO matrix elements [2]. As a consequence, it 
is not possible to find an analytic expression for the amplitudes, which covers the whole kinematic region 
and does not contain functions with implicit imaginary parts. We shall sec below that the most convenient 
representation is obtained by subdividing the kinematical plane for each process according to the cuts into 
four sectors. In each of these sectors, we obtain an analytic expression for the amplitudes with all imaginary 
parts made explicit; joining the different regions together, one obtains a result which is continuous, but not 
analytic, across the cuts. 



3 Helicity amplitudes 

The renormalized amplitude \M) can be written as 

\M) =V^(q;g;q) , (3.1) 

where y M represents the lepton current and denotes the hadron current. With q, we denote the incoming 
quark (momentum —pi) in lepton quark scattering or the outgoing antiquark (momentum +P2) in lepton 
gluon scattering. A convenient method to evaluate the helicity amplitudes is in terms of Weyl-van der 
Waerden spinors, which is described briefly in the appendix of [13] and in detail in [34,35]. 

Using the Weyl-van der Waerden spinor calculus, we can express the lepton current for photon exchange 
in terms of the helicities of the incoming l~(p&) and outgoing f>s). Explicitly, 

V^(l-(-p 5 )-,l-(p 6 )-) = eaj B p 6A p 5B ^§ , V^(l-(-p 5 )+,l-{p 6 )+) = eafp 5AP6B ^ , (3.2) 

Pa Pa 



L lh= Rl hh=- e f^hh- (3-3) 



where for photon exchange, 

It is straightforward to account for charged and neutral current exchange by appropriate replacement of the 
couplings, keeping in mind that in these cases left- and right-handed couplings are no longer identical. The 
lepton current for incident anti-leptons follows from charge conjugation. 
The hadronic current 5 M is related to the fixed helicity currents, S Ab , by 

S„(q+;gX;q-) = R} ih V2aj B S AB (q+;g\; q-), (3.4) 

S fl (q-,g\;q+) = L] ih V2 aj B S AB (q-; g\;q+). (3.5) 

As above, the gauge boson coupling is extracted from S Ab , and charged and neutral current reactions follow 
from appropriate substitutions. 

The currents with the quark helicities flipped follow from parity conservation: 

5 ls('?-;.9A;g+) = {S BA (q+;g{-\);q-))* . (3.6) 
Charge conjugation implies the following relations between currents with different helicities: 

s AB(Q\'9\Q^q) = (- 1 )5 y 4 B (gA ? ;gA;gA,j). (3.7) 
All helicity amplitudes are therefore related to the amplitudes with \ q = + and A, = — . Explicitly, we find 

S- (o+-a+-a-) - a(x v z) V ^ V ° V2B 1 3 (x v z) V ^ V ° V2B 1 -,(x v z) 2lCB^PjA 
S AB {q+,g+,q ) - a{x, y, z) + V, *) MM + ^ * *) M (p 3 p 2 )* 

+ ^ z ) ( p J^ lP2 )* {p ^b + P2AB + Psab) ■ (3-8) 

The other helicity amplitudes are obtained from S AB (q+; g+; q—) by the above parity and charge conjugation 
relations. Current conservation implies the following relation between the four helicity coefficients, 

a(x, y, z) - (3(x, y, z) - j(x, y, z) - S(x, y, z) = 0. (3.9) 

•S12 
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Equation (|3.9p suffices to fix S(x,y, z). 

The three remaining helicity amplitude coefficients a, (3 and 7 are vectors in colour space and have 
perturbative expansions: 



for fi = a, f3, 7 where the dependence on (x, y, z) is implicit. 

The helicity coefficients were extracted from a Feynman diagram calculation in dimensional regulariza- 
tion [36,37], with d = 4 — 2e space-time dimensions, by applying suitable projectors to the amplitudes at the 
required order in perturbation theory. Once these projectors are applied to expose the helicity structure, one 
is left with the task of computing a large number of two-loop integrals. Using integration-by-parts [37, 38] 
and Lorentz-invariance [39] identities, these were reduced [40] (making extensive use of MAPLE [41] and 
FORM [42]) to a small number of so-called master integrals, which were derived in [43]. Insertion of these 
master integrals into the contraction of Feynman amplitudes with the projectors then yields the unrenor- 
malized helicity amplitudes [13]. A partial check was made in [44], where two of the seven colour factors 
contributing to the 7* — > 3 partons two-loop helicity amplitudes were derived by direct evaluation of the 
integrals (i.e. not using a reduction to master integrals), yielding full agreement with [13]. 

The 3j amplitudes obtained in [13] are expressed in terms of two-dimensional harmonic polylogarithms 
(2dHPLs). The 2dHPLs are an extension of the harmonic polylogarithms (HPLs) of [45]. Numerical routines 
providing an evaluation of the HPLs [46,47] and 2dHPLs [48] are available. The implementations apply to 
HPLs of arbitrary real [45, 46] or complex [47] arguments, but only for a limited range of arguments of the 
2dHPLs. This range of arguments corresponds precisely to the 3j kinematics. Formulae for the analytic 
continuation of HPLs and 2dHPLs to ranges of argument relevant to all 2 — > 2 scattering kinematics were 
derived in [49] . Given the restricted range of arguments for the numerical implementation of the 2dHPLs [48] , 
in [49] transformations to map all 2 — > 2 scattering kinematics into this range of arguments were derived. 
These analytical continuations were subsequently confirmed by direct numerical evaluation [50] of the relevant 
master integrals [43] in the relevant kinematical regions using the numerical Mellin-Barnes technique. Using 
the results of [49] , we can perform the analytic continuation of all the unrenormalizcd helicity amplitudes from 
the 7* — * 3 partons kinematics considered in [12] to the kinematical situation of DIS-(2 + l)-jet production. 

Subsequently, renormalization is performed in the MS-scheme by replacing the bare coupling ao with the 
renormalized coupling a s = a s (p 2 ) , evaluated at the renormalization scale 



S e = (47r) e e" £7 with Euler constant 7 = 0.5772 . . . 

and [1q is the mass parameter introduced in dimensional regularization [36, 37] to maintain a dimensionless 
coupling in the bare QCD Lagrangian density; /3q and (3\ are the first two coefficients of the QCD /3-function: 




(3.10) 




(3.11) 



where 



Po 



IIC A - AT R N F 
6 



Pi 



17C\ - 10C A T R N F - QC F T R N F 
6 



(3.12) 



with the QCD colour factors 



C A = N, 



C f 



N 2 



1 




(3.13) 



2N 



For the remainder of this paper we will set the renormalization scale /r 
This procedure yields the renormalized tensor coefficients, 

q(0) = n (0),un 





(3.14) 
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The full scale dependence of the renormalized helicity coefficients can be recovered from the renormalization 
group equation: 



a s (p 2 



2n 



n CD + ( 3^ n (x) + WtfO) )i n (tL)+ 21± n (o) ln 2 ( al 



3/?o 2 



2 / P 



+ 0(o»H. (3.15) 



2tt 

After performing ultraviolet renormalization, the amplitudes still contain singularities, which are of in- 
frared origin and will be analytically cancelled by those occurring in radiative processes of the same order. 
Catani [51] has shown how to organize the infrared pole structure of the one- and two-loop contributions renor- 
malized in the MS scheme in terms of the tree and renormalized one-loop amplitudes. The same procedure 
applies to the tensor coefficients: 



Q« = i«(e)tt (0) +r>« 



finite 



Q< 2 ) = [ —I«(e)I«(e) - ^J«(e) + e' ^)' 2 "} + k] J«(2e) + H^(e) ) to® 
\ 2 e r(l-e) V e / / 



2 e 
+jW(e)12 (1) +0< 2 )' finitc 

where the constant K is 

K = 



67 

18 6 



C A 



10 
~9 



T R N F . 



(3.16) 
(3-17) 



For this particular process, there is only one colour structure present at tree level which, in terms of the gluon 
colour a and the quark and antiquark colours i and j, is simply T° . Adding higher loops does not introduce 
additional colour structures, and the amplitudes are therefore vectors in a one-dimensional space. Similarly, 
the infrared singularity operator I^ 1 ' (e) is a scalar in colour space and is given by 



J (1) (e) = - ; 



2T(1 - e) 
where (since we have set /j, 2 



• , 1 3 fa 
N 1 h — — 

1 e 2 4e 2Ne 



1 ( 1 



(s 13 + s 23 )-- [-, + - s 



iV 



2c 



212 



^123 



(3.18) 



(3.19) 



Using the invariants Sy as defined in Section [2j phase factors distinguishing initial and final state partons 
become obsolete in this expression. Note that on expanding Sy, imaginary parts are generated, the sign of 
which is fixed by the small imaginary part +i0 of all sy . 

Finally, the term of Eq. (|3.16p that involves H i ~ 2 \e) produces only a single pole in e and is given by 



H^ 2 \e 



s«7 



4eT(l - e) 



(3.20) 



where the constant is renormalization-scheme-dependent. As with the single-pole parts of I^'{e), the 
process-dependent can be constructed by counting the number of radiating partons present in the event. 
In our case, in the MS scheme: 



At leading order 



4Ca 



589 
+ 432 ~ 
19 tt 2 
18 + 36 



11tt z 
^72~ 



N 



NNp 





v 

2 C3_ 


41 


-) 


(- 


54 ~ 


48 J 


1 


2 x 

7T 


\ Np 


5 


54 


~ 24 


) N 


+ 27 



-SCa- 



ie 



T 



i 

iV 2 



N 2 . 



a^(x,y,z)=p^(x,y,z) = l 



and 



(3.21) 
(3.22) 
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The renormalized next-to-leading order helicity amplitude coefficients can be straightforwardly obtained to 
all orders in e from the known exact expressions for the one-loop master integrals. For practical purposes, 
they are needed through to 0(e 2 ) in evaluating the infrared-divergent one-loop contribution to the two-loop 
amplitude, while only the finite piece is needed for the one-loop self-interference. They can be decomposed 
according to their colour structure as follows: 

^• Rnn %x,y,z) = Nan(x,y,z) + ^bn(x,y,z) + /3 cn(x,y,z) . (3.23) 

The expansion of the coefficients through to e 2 yields HPLs and 2dHPLs up to weight 4 for an, and 
up to weight 3 for cq. In order to obtain expressions suitable for numerical evaluation (which requires the 
arguments of the 2dHPLs to be in a restricted range) and to make all imaginary parts explicit, we use the 
following decomposition [2] for the coefficients I — a, b, c of the individual colour factors: 

lu(x,y,z) = G(x)Q(l-x)9{y)G(l-y)0{-z)l UA (x + z,-z) 

y + z y 



+e{x-l)e(y)Q(l-y)e(-z)l n ,2 



+e(x)e(i-x)e(y-i)e(-z)iu, 3 



X X 

y + z x 

v ' y 



+9(x - 1)9(2/ - l)e(-z)in l4 ( V±l,— ) (3.24) 

z z 



for the lepton-quark scattering process (|2.ip and 



ln(x,y,z) = G(z)G(l-z)e(y)Q(l-y)G(-x)ln, 5 (x + z,-x) 

x + y y 



+G{z-l)G(y)Q(l-y)e(-x)ln t6 
+6(2)6(1 - z)Q{y - l)6(-z)Zn,7 



z z 
x + z z 



y y 



+0(z ~ 1)9(2/ - l)9(-x)/n,8 ( — , ~- ) (3.25) 

x x ' 



for the lepton-gluon scattering process 

The coefficients ln,i(r, s) are then expressed in terms of HPLs of argument s and 2dHPLs of argument r, 
with s featuring in the index vector. The above choices of argument ensure that < r < 1 — s and < s < 1, 
as required for the numerical evaluation of the 2dHPLs [48]. The explicit expressions are of considerable size, 
such that we shall not quote them here. An example of the size and structure of those coefficients can be 
found in [12], where we explicitly list the helicity-averaged one-loop times one-loop and tree times two-loop 
matrix elements for e + e~ — > 3 jets. It should be noted that these finite pieces of the one-loop coefficients 
can equally well be written in terms of ordinary logarithms and dilogarithms, see [52,53]. The reason for 
expressing them in terms of HPLs and 2dHPLs here is their usage in the infrared counter-term of the two-loop 
coefficients, which cannot be fully expressed in terms of logarithmic and polylogarithmic functions. 

The finite two-loop remainder is obtained by subtracting the predicted infrared structure (expanded 
through to C(e )) from the renormalized helicity coefficient. We further decompose the finite remainder 
according to the colour structure, as follows: 

n (2),fimte (a , ; y) z) = N iA n (x, y, z) + B n (x, y, z) + j^C n (x, y, z) + NN F D n {x, y, z) 

+^E n (x,y,z) + N 2 F F n (x,y,z) + N Fy (Jl - N^j G n {x,y,z) , (3.26) 

where the last term is generated by graphs where the virtual gauge boson does not couple directly to the 
final-state quarks. This contribution is denoted by Np.y and is proportional to the charge weighted sum of 



G 



the quark flavours. In the case of purely electromagnetic interactions we find, 

N F>1 = . (3.27) 

Including electroweak vector boson interactions, the same class of diagrams yields not only a contribution from 
the vector component (obtained by appropriate replacement of the couplings in the above formula [13]), but 
also a contribution involving the axial couplings of the vector boson under consideration. This contribution 
vanishes however if summed over isospin doublets. 

We apply the same decomposition as used at one loop for the coefficients L = A, . . . , G of the individual 
colour factors: 

L Q (x,y,z) = G(x)G(l-x)G(y)G(l-y)Q{-z)L n ,i(x + z,-z) 

y + z y 



-Q(x-l)Q(y)Q(l-y)Q(-z)L n ,2 



+0(1)6(1 - x)0(y - l)G(-z)L n , 3 



X X 

y + z x 

y ' y 



+9(x - 1)9(2/ - l)0(-*)Ln,4 ~\ ( 3 - 28 ) 

for the lepton-quark scattering process (|2.ip and 

L n (x,y,z) = 0(z)0(l - z)G(y)G(l - y)G(-x)L n , 5 {x + z, -x) 

x + y y 



+Q(z - l)0(i/)0(l - y)G{-x)L nfi 
+0(2)0(1 - z)G{y - l)G(-x)L nj 



z z 
x + z z 



y y 



+G(z - 1)0(1/ - l)G(-x)L n>8 ( Z±l 1 ) (3.29) 

x x ' 



for the lepton-gluon scattering process 

The helicity coefficients contain HPLs and 2dHPLs up to weight 4 in the A, B, C, G-terms, up to weight 
3 in the D, _E-terms and up to weight 2 in the F-term. The size of each helicity coefficient is comparable to 
the size of the helicity-averaged tree times two-loop matrix element quoted in [12]. We therefore refrain from 
quoting them explicitly. The complete set of one-loop and two-loop coefficients in FORM and FORTRAN 



format is provided with the sources of the paper on the archive |http://arxiv.org 



4 Conclusions and Outlook 

In this paper, we have applied the analytic continuation procedures for two-loop four-point functions with one 
off-shell leg derived in [49] to obtain the two-loop helicity amplitudes relevant to DIS-(2+l) jet production 
from the amplitudes for 7* — > qqg derived earlier. Besides enabling the computation of the virtual two-loop 
corrections to the DIS-(2+l) jet cross section and related event shape observables (which follow by interfering 
the amplitudes with their tree level counterparts and summing over all external helicities), knowledge of the 
helicity amplitudes provides additional information about the scattering process. In particular, oriented event 
shapes as well as jet rates in polarized DIS can be computed from them. These matrix elements represent one 
of the missing ingredients for NNLO computations of deep inelastic jet production. Such calculations require 
the construction of a parton-level event generator program, which contains all partonic channels contributing 
to the observable under consideration. Since the one-loop and tree-level contributions at NNLO contain more 
partons than hard jets, they will contribute infrared real radiation singularities, which have to be subtracted in 
the corresponding partonic channels, and added to the two- loop virtual contribution to yield finite expressions 
suitable for numerical evaluation. While a formalism for NNLO subtraction is well established for processes 
involving only final-state partons, some extra work is required to subtract singularities arising from unresolved 
emission off initial-state partons. Work in this direction is ongoing. 
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